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Abstract 

We investigate the renormalizability of the classical 4 theory at 
finite temperature. We calculate the time-dependent two point func- 
tion to two loop order and show that it can be rendered finite by the 
counterterms of the classical static theory. As an application the clas- 
sical plasmon damping rate is found to be 7 = A 2 T 2 /15367rm. When 
we use the high temperature expression for m given by dimensional 
reduction, the rate is found to agree with the quantum mechanical 
result. 
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1 Introduction. 



In recent years there has been an increasing interest in real-time phenomena in 
finite temperature quantum field theory. The complexity of the phenomena of- 
ten calls for numerical simulations. Because it is extremely difficult to simulate 
quantum mechanical real-time quantities, approximations have to be made. In 
particular the classical approximation has been used to simulate diffusion of the 
Chern-Simons number in gauge theories [1-4] and the dynamical properties of the 
electroweak phase transition ||. It has been argued that the classical descrip- 
tion can be accurate for physical quantities which are mainly determined by low 
momentum modes of the theory, because these modes will be highly occupied at 
high temperatures. However, a clear demonstration is complicated by the fact 
that there may be modes with momentum higher that the temperature. Exclud- 
ing these by a cutoff in an effective field theory allows for a controlled derivation, 
but it leads to very complicated dynamics with nonlocal noise terms |§. It is 
much simpler to use cutoffs much higher than the temperature and deal with 
dynamics which is local at this cutoff scale. In this letter we shall study the limit 
of infinite cutoff of the (noiseless) classical theory in perturbation theory, taking 
the (/> 4 theory as an example. 

Classical time dependent correlation functions are determined by solving clas- 
sical equations of motion for arbitrary initial conditions, and performing a Boltz- 
mann weighted sum over these initial conditions. For time independent correla- 
tion functions, this prescription reduces to calculations in a three dimensional, 
euclidean field theory. In fact, this static theory is of the same form as the di- 
mensional reduction approximation in finite temperature field theory [0, ||. In 
one method of dimensional reduction, the values of the coupling constants are 
precisely specified by a matching procedure between the four and the three di- 
mensional theory [^]. In the classical theory the coupling constants are apriori 
arbitrary. The dimensionally reduced theory has been used extensively for the 
study of the electroweak phase transition 0. 

The euclidean theory coming from dimensional reduction is superrenormal- 
izable and therefore, time independent correlation functions can be made finite 
by merely a mass renormalization. It is not apriori clear that time dependent 
correlation functions can be made finite by the same renormalization. In |J it is 
shown for a 4 -theory that the time dependent two point function can be made 
one loop finite by supplying the mass counterterm of the static theory. However, 
at two loops the situation is different because then the correlation functions of 
the static theory mix with time dependent Green functions. We shall study this 
question by calculating the classical time dependent two point function to two 
loop order. 

We will solve the equations of motion by applying naive perturbation theory, 
which leads to resonant or secular terms in each step of the expansion. These 
can be avoided by adding counterterms to eliminate the secular terms. We will 
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show that the counterterms needed to make the theory finite can be chosen such 
that they also partially take care of the resonances. 

We also calculate the classical plasmon rate and show that it is equal to the 
quantum mechanical expression at high T, when the connection is made with 
dimensional reduction. 

2 Definition of the correlation functions. 

We consider a scalar field theory at temperature T = 1/(3, with the following 
hamiltonian 

E = J** fa* + ^(V0) 2 + + ^ + e) . 

Classical real time n point functions are defined in the following way 

(0( Xl , h)... 0(x n , t n )) = | J VnV^e-^^Uixu h) . . . 0(x n , t n ), 

where 

Z = J VnVfa-W*'® , 
and 0(x, t) is the solution of the equations of motion 

0(x,t) = 7r(x,t), 

7r(x,t) = V 2 0(x,t) 
with initial conditions 

0(X, t ) = 0(x), 7r(x, t ) = 7T(x). 

The integrals over 7r and are performed at t = t (summing over the initial 
conditions). Since we consider equilibrium correlation functions, t is arbitrary. 
In particular the two point function will only depend on ti—t 2 . Choosing infinite 
volume it depends also only on x x — x 2 . In the following we choose parameters 
such that the theory is in its symmetric phase. 

For equal times t± — . . . — t n — t the 7r-field decouples and we are left with 
a three dimensional theory for <fi. This three dimensional theory is superrenor- 
malizable and it can be made finite with merely a mass renormalization. The 
counterterms are determined by the leaf diagram and the momentum independent 
part of the setting sun diagram. Writing 

fx 2 = m 2 — 5m 2 , 5m 2 = \m\ + \ 2 m\, 



-/i 2 0(x,t)--0 3 (x,t), 
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a) b) c) 

Figure 1: Propagators, a) So(x — x', t — t'), b) So(x — x',0), c) G^(x — 



*,t-1f). 



we will use 



ml 



rrir, = 



-T 



d 3 k 1 
(2^ 



At AT - — mT, 
4tt 2 8tt 



(1) 



-T 



d 3 k x d 3 k 2 d 3 k 3 (2vr) 3 (5(k 1 + k 2 + k 3 ) 
J^nf J2nf J2nf ~ 



T 2 A 2 
32^ log(-)+ finite, 



(2) 



with A the momentum cutoff and cu^ = k 2 + m 2 . The choice of the finite parts 
of the counterterms implied by the above expressions will turn out to be conve- 
nient for the moment. With this mass renormalization the equal time correlation 
functions are finite. 

The first few terms of a perturbative solution = O + A0i + A 2 02 + • • • to 
the equations of motion can easily be constructed. The zeroth order solution is 
given by 



0o(x,£) 



d 3 k 



,ik-x 



7r(k) 

0(k) cos(cj k (t - t )) H sin(cu k (t - t )) 



where 0(k) and 7r(k) are the Fourier components of the initial conditions. Inte- 
gration over the initial conditions gives the free two point function 

(0 o (x,t)0 o (x',t'))o-^o(x-x',t-t')=T / ^Hx-xO eosK(t-^ 

Imposing the initial conditions <pi = 0i = 02 — <p2 — • • • — the higher order 
terms can be found with the help of the retarded Green function 



G£(x - x', t - t') 



d 3 k eik . (x _ x0 ^^ _ t ,^sm(u k (t - t')) 



(27T) 



The solution to second order is given by 

0(x,t) = O (x,t) 

- A J dt' J d 3 x' G*(x - x', t - t') 
+ A 2 J dt' J d 3 x' G^(x - x', t - t')mlMx, t 



-^(x',O-m 2 0o(x',O 



Figure 2: First order contribution to the two point function, the first order coun- 
terterm is denoted with a cross. 



+ A 2 J dt' J d 3 x' J dt" J d 3 x" G*(x - x', t - t')G*{x' - x", t' - t") 



2! 



^(x',0 



<Pl(x",t")-m(Mx",t") 



3! 



where the integrals over t',t" start at to- We expressed fi 2 in the equation of 
motion in terms of m 2 and 5m? and expanded in A, using the same coefficients m\ 
and m| as found earlier for the equal time case. Expanding also the Boltzmann 
weight, all correlation functions can be expressed in terms of the free n-point 
functions of O ( X , t), which are just products of the two point function S Q , by 
Wick's theorem. The various propagators are given in figure |I[ It is not clear at 
this point that the static counterterms suffice also for making the time dependent 
correlation functions finite, but we shall see in the following that this is indeed 
the case. 



3 First and second order calculation of the two 
point function. 

Using the framework set up in the previous paragraph, we calculate the time 
dependent two point function to first order in A. The result is 

(0(x 1; ti)0(x 2 , t 2 )) = So(Xi - X 2 , ti - t 2 ) 

- |/3A j rf 3 x5 (x 1 -x,t 1 -t )(5o(0,0)-2m 2 )5 (x-x 2 ,t -t 2 ) 

-^A J dt y d 3 xG^(xi -x,ti -t) (£0(0,0) -2m 2 )S , (x-x 2 ,t-t 2 ) 
+ (xi,ti) <-> (x 2 ,t 2 ) 

The corresponding diagrams are given in figure 0. Performing the time integral, 
the correlation function becomes in momentum space 




T cos(a;pt) 
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AT / [ d 3 k 1 ,\ 1 . . . . ,. 

~4~\J (2^)3^f ~ 2m 7 M sm K*) + 2cos(w p t)] 



where t = t\ — t 2 . We see that the one loop divergence in this time dependent 
correlation function is cancelled by the one loop counterterm (j^) of the static 
theory. The particular choice (|l|) for the finite terms also cancels the resonance 
due to the secular term. 

We now turn to the second order calculation of the correlation function. Using 
the result for the first order counterterm, which cancels completely all the one 
loop contributions, the contribution to second order is 

(</>(xi, £i)0(x 2 , t 2 )) = So(xi - x 2 , £i - t 2 ) 

(a) +^P 2 \ 2 J d 3 xd 3 x' 5 (x! -x,t! -£ )^(x-x , ,0)5 (x / -x 2 ,t -t 2 ) 

(b) + (3m\\ 2 J d 3 x S (xi - x, t\ — t )S {x - x 2 , t - t 2 ) 

(c) + i/?A 2 J dt J d 3 x d 3 x' G?J(xi - x, tj. - t)5g(x - x', t - t„) 

xS'o(x / - x 2 ,t - h) + (xi,ti) <-> (x 2 ,t 2 ) 

(d) + ^A 2 y dt dt' y d 3 x d 3 x G^(xx - x, ti - t)G?(x - x', t - t') 

xS 2 (x-x',t-£ , )^o(x / -x 2 ,t'-t 2 ) + (x 1 ,t 1 ) «-> (x 2 ,t 2 ) 

(e) +m 2 A 2 y y d 3 xG^(xi -x,ti -£)S , (x-x 2 ,£-t 2 ) 

+ (xi,ti) <-> (x 2 ,t 2 ) 

(f ) + ^A 2 y dt' y d 3 x d 3 x' G^( Xl - x, *i - t)S 3 (x - x', t - 

xG*(x 2 -x',£ 2 -t'). 

The diagrams are given in figure |3|. Contributions (a) and (b) come from the 
second order expansion of the Boltzmannfactor, (c) from a combination of the 
first order expansion of the Boltzmannfactor and the first order solution of the 
equation of motion, (d) and (e) from the second order solution and (f) from 
the combination of two first order solutions to the equations of motion; (a) and 
(b) are the normal equal time contributions, (c), (d), (e) and (f) vanish when 
t\ = t 2 = to- 



4 Finiteness of the second order contribution. 

We will now see if the counterterm can be chosen in such a way that the second 
order contribution is finite. First we go over to three dimensional momentum 
space. Approximating in expressions (a) to (f) the sines and the cosines by 1, we 
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a) 



b) 



c) 





d) 



e) 



f) 



Figure 3: Second order contribution to the two point function; here the cross 
denotes the second order counterterm. 



see that there are possible linear ((d)) and logarithmic ((a), (c), (f)) divergences 
(actually, the indication of a linear divergence disappears after doing the time 
integrals). Inspection then shows that the momentum dependent part is finite, 
such that only the momentum independent part needs to be studied. We will take 
the external momentum p to be zero, the internal momenta are labelled as ki, k 2 
and k 3 = — k x — k 2 . Using the fact that the correlation function only depends on 
ti — t 2 , we take £ 2 — to — and t\ = t. This immediately sets to zero half of the 
contribution from (c), (d), (e) and the complete contribution from (f). Doing the 
time integrals, the expressions ± c<j kl ± c<j k2 (with all four combinations of the 
±'s) often pop up, and we will denote these by 

^ = ei^kj + e2^k 2 + e 3 ^k 3 , 

with €i — ±1. 

Performing the time integrals, using the symmetry of the resulting expres- 
sion to replace SQe^u^ — > Q 2 and combining several contributions, we find the 
following result, 



S(0,t) 
(a+b) 

(d+e) 

(c+d) 



Tcos(mt) 



+ 



+ 



+ 



TX 2 cos(mt) 

m 1 

TX 2 tsm(mt) 
2m~ 3 

t 3 a 2 r d 3 h 



II 

6 

-rp2 
IT 



d 3 ki d 3 k 2 d 3 k 3 (2n) 3 5(k 1 



k 3 ) 



(2tt) 3 (2tt) 3 (2tt) 3 
d 3 h d 3 k 2 d 3 k 3 (27r) 3 5(k 1 + k 2 + k 3 ) 



(2tt) 3 (2tt) 3 (2tt) 3 
d 3 k 3 (2ir) 3 5(k l + k 2 + k 3 ) 



2 2 2 



m 



+ m 2 



48 J (2tt) 3 (2tt) 3 (2tt) 



^ kl ^k 2 ^k 3 



x 



^ Q 2 

ei£2£3 



t sin(mt) cos(fii) — cos(mt) 



2m 



+ 



n 2 



(3) 
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In parenthesis we have indicated which diagrams contributed. Both (a+b) and 
(d+e) give the same result for the second order counterterm, it is given by @. 
This is an important result because it shows that the counterterms of the static 
theory are sufficient to make the time dependent two point function two loop 
finite. With the choice (fj) (again in particular the finite part), the above (a+b) 
and (d+e) are zero. The contribution indicated by (c+d) is ultraviolet finite. 
The possible divergences for Q — > ±m are actually not present. For example, for 
n — > m, 



1 



Q2 _ m 2 



tsin(mi) cos(fit) — cos(mt) 



2m 



n 2 — m 2 



t sin mt t 2 cos mt 



8m 3 



8m 2 



+ 0{Q 



m) 



Let us record here also the fact that the two loop correction is 0(t 4 ) for t — > 0. 



5 Damping rate 

As an application we derive an expression for the damping rate 7, assuming that 
the dominant large time behavior of the correlation function at zero momentum 
is given by Z(T/m 2 ) exp(— jt) cos((m + 8)t). Since 7 and Z — 1 are 0(A 2 ) and 
5 at least 0(A) we can identify 7 from the coefficient of — tcos(mt) in the large 
time behavior of ([3[). In terms of 



A 2 
6 



e /n 



eie2<:3 " 3=1 
\4 



T 



X2vr) 3 2u; kj w k . 
(27r) 4 5(ki + k 2 + k 3 )5(e 1 u kl + e 2 ^ k2 + e 3 ^k. 
we can rewrite the 0(X 2 ) contribution to (|3|) in the form 

1 



S {2) (0,t) 



m A 



t sin mt cos Vtt — cos mt 
+ 



2m 



9? 



The function w(Q) is even in Q — > — Q, continuous and vanishes like at large 
Q. A large time behavior cx t cos mt can only come from the singularities at 
Q = ±m, any other portion of the Q integral gives subdominant behavior. Hence 
we may conveniently replace w(Q) —>■ w(m) and evaluate the resulting integral, 



S (2) (0,t) 



from which we identify 



w(m\ 



7 



sin mt t cos mt 



Am 3 



w(m) 



Am 2 



AT 



The e's giving nonvanishing contributions to w(Q) are ei = e 2 = 1, €3 = — 1 
(cycl) for Q > and similar for Q < 0. Hence w(Q) has the form of a phase 



8 



space integral for the scattering (0, ft) + (k^cc^) — > (ki^kj + (k2,^k 2 ) etc. 
The particles are distributed according to T/u, which is the high temperature 
form of the quantal distribution [exp(u/T) — 1] . Evaluation of w(m) gives for 
the damping rate 

7 = ■ (4) 

15367rm 

6 Matching to the quantum theory 

The classical theory can be seen as a high temperature approximation to the 
quantum theory. Then the parameters in the classical theory are effective pa- 
rameters which are to be chosen such that the classical correlation functions 
match those in the quantum theory as well as possible. For the static theory the 
approximation is equivalent to a form of dimensional reduction and this makes it 
possible to determine the parameters in the potential energy part of the hamil- 
tonian in terms of the parameters in the quantum theory M. General matching 
rules are obtained in ref. |7|. For example, matching the cosmological constant, 
mass and coupling parameters to one loop order we find 

tt 2 T 4 m 2 T 2 m 4 , v 2 A 3 T / A 1\ ATm 2 
+ ~ Trr^T m — - I In — - - I r^-, 



90 24 64vr 2 v 2 6vr 2 V T 3/ 4vr 2 

, 2 T (T 2 m 2 , z/ 2 AT 

u = m + A In — =- 

P V24 32vr 2 v 2 4tt 2 

A = A ( 1 - In 

where vt = 47rexp(— je)T and less dominant terms as T — > oo are suppressed. 
The quantities on the right hand side are dimensionally renormalized parameters 
in the MS-bar scheme and v is the corresponding mass scale. Notice how the 
Rayleigh- Jeans divergence is canceled by the A 3 counterterm, leaving the familiar 
quantum expression for the pressure as just a finite renormalization in the classical 
theory. 

The above relation for fi 2 now implies with (|I|), choosing the MS-bar scale 
v = u T , 

2 2 JT 2 fflT 
m = m + A — - 



24 8tt 

Inserting this relation in the result (J|) for the damping rate gives for large T 

A 3/2 T 

7 



128^6 tt' 

which is indeed the plasmon rate obtained in the quantum theory [|ITJ. 
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7 Discussion 



Among the time-independent correlation functions only the two-point function is 
superficially divergent, because of the superrenormalizability. Inspection of the 
time-dependent correlation functions leads to the conclusion that, after perform- 
ing the intermediary time integrals and replacing the remaining sines and cosines 
by 1, the overall degree of divergence of the n > 2-point functions is less than zero. 
Only the n = 2 case is then non-trivial and we have shown to second order in 
A that the time-dependent classical two-point function is finite, after taking into 
account the counterterms of the static theory. We conclude therefore (but have 
of course not given a proof) that with these counterterms all finite temperature 
time dependent correlation functions in 4 theory are cutoff independent. 

The coupling A was taken to be cutoff independent, consistent with the super- 
renormalizability of the static theory. We would like to remark that this static 
three dimensional theory is quite non-generic from the Statistical Physics point 
of view. The coupling A3 of a three dimensional theory has the dimension of 
[mass] and the dimensionless coupling at scale v is \(y) = A 3 (z/)/V. We may 
call \b = A (A) the bare coupling and Xr = A(m) the renormalized coupling. 
Since the system gets critical as the ratio A/m approaches infinity, one might 
expect that the renormalized coupling approaches the nontrivial infrared stable 
fixed point A* of Ising-like systems (the Wilson-Fischer point). This is not the 
case here, because Xb = AT/A is non-generic: it is is automatically tuned to 
the gaussian fixed point as A/m —>■ 00. The renormalized coupling \r « XT/m, 
which can be arbitrary in the interval (0, A*). 

For a clearer understanding of the classical theory as an approximation to the 
quantum theory one would like to go beyond the usual dimensional reduction 
involving only time independent quantities. Matching time dependent quantities 
may introduce at two loop order also a finite renormalization of the coefficient 
of 7r 2 in the hamiltonian. We have not yet carried out such time dependent 
matching. For a different approach to time dependent quantities, based on an 



expansion in Planck's constant h, see ref. [Il . 

The fact that the analytic expression for the classical plasmon damping rate 
combined with parameter matching according to dimensional reduction leads to 
the correct quantum mechanical result is quite encouraging for the physically 



more interesting case of the SU(2)-Higgs model (cf. ref. [12] where eq. 4.3 gives 
the impression that the classical gluon damping rate is uv-divergent; quantumme- 
chanically the gluon damping rate is of course uv-finite but subtle in the infrared, 



see e.g. |13[). Numerical simulations in the classical SU(2)-Higgs model clearly 
show the damping phenomenon and appear to be compatible with lattice spacing 
independence of the rates |Tj|. 
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